Reduced effective spin-orbital degeneracy and spin-orbital ordering in paramagnetic 

transition metal oxides: Sr 2 Ir04 vs. Sr 2 Rh04 
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We discuss the notions of spin-orbital polarization and ordering in paramagnetic materials, and 
address their consequences in transition metal oxides. Extending the combined density functional 
and dynamical mean field theory scheme to the case of materials with large spin-orbit interactions, 
we investigate the electronic excitations of the paramagnetic phases of SralrCU and Sr2Rh04. We 
show that the interplay of spin-orbit interactions, structural distortions and Coulomb interactions 
suppresses spin-orbital fluctuations. As a result, the room temperature phase of Sr2lr04 is a para- 
magnetic spin-orbitally ordered Mott insulator. In Sr2Rh04, the effective spin-orbital degeneracy is 
reduced, but the material remains metallic, due to both, smaller spin-orbit and smaller Coulomb in- 
teractions. We find excellent agreement of our ab-initio calculations for Sr2Rh04 with angle-resolved 
photoemission, and make predictions for spectra of the paramagnetic phase of Sr2lr04. 
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Probably the most important consequence of relativis- 
tic quantum mechanics in solids is the coupling of spin 
and orbital degrees of freedom. The concept of spin-orbit 
coupling (SOC) has been known for more than half a cen- 
tury, and its importance in magnetic materials has been 
recognized early on. Nevertheless, SOC in non-magnetic 
materials has only recently become a hot topic in con- 
densed matter physics, in the context of topological in- 
sulators [l|, within the search for systems for quantum 
computing applications. Still, the interplay of spin-orbit 
(SO) interactions with electronic Coulomb interactions 
in paramagnetic materials remains a largely unexplored 
field. This is at least partially due to the seemingly mutu- 
ally exclusive regimes of their acting: SO interactions are 
strongest in heavy elements, that is in transition metal 
compounds with 5d (and, to a lesser extent Ad) electrons. 
Coulomb interactions, on the other hand, are believed to 
be most efficient in 3d compounds, due to the more local- 
ized 3d orbitals. It is only recently that this separation 
has been questioned, e.g. in [2-5]. 

In this Letter, we present a particularly striking ex- 
ample for the interplay of SOC and Coulomb interac- 
tions: in the room-temperature phase of the 5d transi- 
tion metal oxide (TMO) Sr2lr04 the combined effect of 
SOC and distortions from the ideal K2NiF4-structure is 
strong enough to suppress the effective degeneracy of the 
coupled spin-orbital degrees of freedom. Even moder- 
ate Coulomb interactions then induce a Mott localized 
state. The resulting state - despite of being paramag- 
netic and paraorbital, i.e. in the absence of any mag- 
netic or orbital order - displays a "spin-orbital order" 
in the sense that the only hole in the t2 ff manifold oc- 
cupies a state of well-defined t2 9 -projected total angular 
momentum J 2 . This notion generalizes the concept of 
orbital order, well-known in TMOs, to the case where 
neither spin nor orbitals are good quantum numbers. We 



present - entirely from first-principle calculations - a sce- 
nario for the paramagnetic insulating state of Sr2lr04, 
in comparison to the isostructural and isoelectronic 4d 
compound Sr 2 Rh04. For the latter, we find that the 
smaller SOC and Coulomb interactions induce a partially 
spin-orbitally polarized metallic state. We calculate pho- 
toemission spectra for both materials and find excellent 
agreement with available experiments. 

Sr2lr04 is a 5<i-TMO with a tetragonal crystal struc- 
ture whose symmetry is lowered from the K2NiF4-type, 
well-known in Sr 2 Ru04 or La 2 Cu04, by 11° rotation 
of its IrO@ octahedra around the c-axis as in its 4d- 
counterpart Sr 2 Rh04 [a, [Z(. Although each Ir site acco- 
modates 5 electrons, Sr2lr04 exhibits an insulating be- 
havior at all temperatures with an optical gap of about 
0.26 eV at room temperature [8(- Below 240 K, canted- 
antifcrromagnetic (AF) order sets in, with an effective 
paramagnetic moment of 0.5/is/Ir and a saturation mo- 
ment of 0.14ur/Ix 9jl. This phase has triggered much 

Y\ I — I I — 1 1 — 
experimental and theoretical work recently [3j, |4j, I lOt lll| , 

highlighting in particular the importance of the SOC Cso- 
SOC was also shown to be relevant in Sr 2 Rh04, which is 
a paramagnetic metal down to 36 mK [12j : density func- 
tional calculations within the local density approximation 
(LDA), augmented by Coulomb interactions within the 
LDA+U method, reproduce the Fermi surface only if SO 
interactions are taken into account 
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Here, we focus on the paramagnetic insulating phase 
of Sr2lr04 above 240 K, which has not been addressed 
by theory before. We analyze the electronic properties 
by a combination of the LDA with dynamical mean field 
theory (DMFT). Our method is a generalization of the 
LDA+DMFT scheme as implemented in [15j based on 
the Wien2k package jiff , extended to include SO interac- 
tions [17| . The Coulomb interactions are calculated from 
the constrained random phase approximation (cRPA) 
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Figure 2: LDA band structure of Sr2ir04 (left) and Sr2Rh04 
(right), projected on the j e ft=l/2 (top), j c a=3/2 \rrij\=3/2 
(middle), and j c ff=3/2 \rrij\ — 1/2 (bottom) spin-orbitals. 



Figure 1: LDA band structures of Sr2lr04 without distortions 
and without SOC (a), with distortions without SOC (b), with- 
out distortions with SOC (c) and with both, distortions and 
SOC (d). 



18], taking matrix elements in the same set of local- 
ized orbitals as used in the DMFT calculations [19j. For 
Sr2lrC>4 we estimate the Slater integrals as F°=2.2eV 
and J=0.3eV, giving an intraorbital Hubbard U in the 
half-filled jeff=l/2 orbitals of f/=2.25eV 20]. 

Within an LDA band picture, Fig. Q] (d), a metallic 
solution - at variance with experiments - is obtained for 
Sr 2 Ir04. We construct Wannier functions for the t 2g or- 
bitals from the LDA bands within the energy window 
[-3.0, 0.5] eV. The large SOC of about Cso«0.4eV splits 
these states into a quartet of states, commonly labeled as 
j e ff=3/2, and a higher lying doublet j e s=l/2. Each state 
is twice degenerate in ±toj. The four j c g=3/2 states are 
almost filled; we find ?t-3/2,|i/2|=1-98 and %/2,|3/2|=l-84 
for the j e ff=3/2 orbitals. The j e ff=l/2 states thus slightly 
exceed half-filling with n 1 / 2 =l.l6- In the left panel of 
Fig. [21 we plot the j e s character of the Wannier orbitals. 
The four bands that cross the Fermi level are not purely 
formed by the j e ff=l/2 orbitals, but there is a slight over- 
lap of the j c ff=l/2 and j e ff=3/2 |mj|=3/2 characters, 
particularly obvious at the T point. Similar conclusions 
were drawn for the AF phase in Ref . [ll| . 

We now redefine the j c ff=l/2 and j c ff=3/2 |mj|=l/2 
atomic states, yielding a diagonal density matrix of the 
local problem. Thus, our definition takes into account 
the structural distortions and a tetragonal crystal field 



in each IrOg octahedra. Our coefficients for the j e e=l/2 
state are similar to those obtained for the AF phase in 
Ref. 0. 

We now turn to our LDA+DMFT results. In Fig. El 
we display the corresponding spectral functions. Within 
LDA+DMFT, supplemented by the cRPA interactions, 
an insulating solution with a Mott gap of the size of 
the optical gap measured at room temperature (about 
0.26 eV [i|) is obtained. Despite the fact that the ARPES 
measurements of Ref. [3j were performed in the AF phase, 
a comparison of the total spectral function to the ex- 
perimental energy distribution curves shows qualitative 
and even quantitative agreement along the direction T-X 
and T-M (between -0.5 eV and -1.5 eV). Along the M-X 
direction a band around -leV is identified. From the 
orbitally- resolved spectral functions (Fig.[3]b and c), one 
can locate the lower j ff=l/2 Hubbard band at about 
-0.5 eV, in agreement with the ARPES data. 

The spin-orbital polarization is enhanced when taking 
into account Coulomb correlations: the Wannier orbital 
jcff=l/2 is now exactly half- filled, and the upper Hub- 
bard band is of j c ff=l/2 type only. We are thus dealing 
with a state exhibiting a "spin-orbital order" in the sense 
of a well-defined j e g quantum number. In contrast to the 
AF phase, where it is also the j ff=l/2 state that carries 
the hole, we do not have an ordering of the correspond- 
ing m., cff =±l/2 quantum number. Also, the orbital oc- 
cupations of the original t 2g orbitals to the j c ff=l/2 spin- 
orbital are roughly equal. We are thus facing the remark- 
able situation of "spin-orbital order" , with neither orbital 
nor magnetic order. 

As it is well-known from studies of multi-orbital Hub- 
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Figure 3: Momentum-resolved spectral functions A(k, uS) of 
the paramagnetic phase of Sr2lr04 from LDA+DMFT at 
T=300K. 



bard models the critical interaction for the formation 
of the insulating state is lowest for the one-band case, 
and increases with degeneracy (for a review see 21|). In 
Sr 2 Ir04, the Hubbard interactions calculated from cRPA 
are large enough to induce a Mott insulating state in 
a half-filled one-orbital but not in a quarter-filled two- 
orbital or a 1/6- filled three-orbital system. We can thus 
conclude that the reduced effective spin-orbital degener- 
acy is the reason for the insulating nature of Sr2lr04. An 
analogous suppression of orbital degeneracy, albeit purely 
induced by cristal field splittings, has been studied in the 
model context in 



22, 23]. and has been found to make 



LaTi0 3 and YTi0 3 Mott-insulating [24 1. 

In Sr2lrC>4 the suppression of spin-orbital fluctuations 
is a consequence of the combined effect of SOC and the 
tetragonal distortions. To demonstrate this, we do a nu- 
merical experiment in which we switch on or off the two 
effects separately and monitor the critical value of U for 
the Mott metal-insulator transition. 

In Fig. Q] (a) both the SOC and the structural dis- 
tortions arc omitted (we use the lattice parameters of 
Rcf. 25]). The band structure of Sr2lr04 is then very 
similar to the one of Sr2Ru04, except for the increased 
filling compared to the d 4 Ru-compound. The t2 ff orbitals 
almost equally accomodate the 5 electrons (n^ =1.612, 
n d x ,,rf ,=1-696). The insulating Mott state is reached 
for F°=3.6eV and J=0.3eV with enhanced orbital po- 



larization. The d xy band gets filled with increasing cor- 
relations, approaching a completely filled d xy band and 
two 3/4-filled d xz / yz bands (which is the atomic ground 
state). 

In Fig. [2(b), the SOC is neglected but the structural 
distortions are included. This lowers the symmetry of 
Si'2lr04 and results in a four-times larger unit cell and 
four-folded bands. Similarly as in Sr 2 Rh04 [l3|, a d xy - 
d x 2_ y 2 hybridization gap opens between the e g and the 
t2 9 bands (between 0.4 and 1.3eV), and the d xy band 
becomes almost filled. The orbital polarization is thus 
enhanced compared to the undistorted case, and the 
Mott transition occurs for smaller values of F°. Indeed, 
LDA+DMFT gives the insulating state for F°=3.0eV 
with J=0.3eV. 

In Fig. Q] (c), we show the band structure of "undis- 
torted" Sr2lr04 but with SOC. Three j e g bands can 
be identified: the j g=l/2 one with an occupation of 
714/2=1-20 lies above the two j c ff=3/2 ones, with fill- 
ings n 3 /2.|3/2|=l-92 and "-3/2,|i/2|=l-96. We find this 
three band system to be insulating at -F°=3.0eV with 
J=0.3eV. Although the problem is effectively reduced 
to an almost half-filled one-band model, the large band- 
width of the j cf f =1/2 band prevents a smaller value for 
the critical F°. 

Finally, panel Fig. [B(d) depicts the band structure of 
the true compound. In this case, the Mott transition oc- 
curs between F°=2.1 and 2.2 eV for J = 0.3 eV, the val- 
ues that have been estimated by our cRPA calculations. 
The combination of both the structural distortions and 
the SOC is thus necessary for Sr2lr04 to be insulating. 

In particular, we would expect Sr2lr04 to be metallic, 
if either the distortions or the SOC were smaller. Such 
a situation may be realized if the material is strained, or 
grown on a substrate [26j. It also occurs in Sr 2 Rh04, the 
isostructural and isoelectronic 4d counterpart of Sr 2 Ir04. 
As expected for 4<i-orbitals, the bandwidth is smaller. 
As before, we construct Wannier j e g orbitals from the 
bands within the energy window [— 2. 67, 0.37] eV; the 
Wannier j c g character along the bands of Sr2Rh04 is 
shown on the right hand side of Fig. [5] Again, there 
is a mixture of the j e ff=l/2 and j e ff=3/2 |mj|=3/2 or- 
bital character along the four bands crossing the Fermi 
level. However, since the SOC is weaker in this com- 
pound (Cso=0.161eV) than in its 5ci-counterpart, the ef- 
fective splitting between the bands j e g=l/2 and j e ff=3/2 
\rrij |=3/2 is smaller and these two states have similar 
filling (m/ 2 =1.42, n 3 /2,|3/2|=l-64), whereas the j ff=3/2 
|toj |=l/2 state is filled (ri3/2,|i/2|=l-96). The system is 
thus close to a 3/4-filled two-band model, which favours 
a metallic state. 

Moreover, the Hubbard interactions are smaller than 
in Sr 2 Ir04: a cRPA calculation as above gives F°=1.6eV 
and J=0.3cV for Sr2Rh04. Indeed, the weaker hy- 
bridization of the Rh-4d states with the 0-2p locates the 
latter about 1 eV higher in energy than in Sr2lr04. As a 
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Figure 4: Momentum-resolved spectral function A(k.,w) 
of Sr2Rh04 (a). Bottom pannels display the orbitally- 
resolved spectral densities for j e ff=3/2 \rrij\=l/2 (b), j e ff=3/2 
\m,j\=3/2 (c) and j eB =l/2 (d) at T=300K. 



result, the Coulomb interactions are screened more effi- 
ciently in Sr2Rh04 compared to Sr2lr04. By computing 
the interaction within cRPA for an artificial Sr2Rh04 
system in which we shift the O-p-states down by 1 eV we 
verify that F° is indeed increased. 

Our LDA+DMFT calculations, Fig. IU are in good 
agreement with ener gy distribution curves obtained by 
ARPES at 10 K J27|, T^l- We oberve a large electron- 
like pocket of radius 0.65 — 0.69 A -1 and a smaller hole 
pocket of radius 0.19 — 0.18 A -1 in the undistorted Bril- 
louin zone. Further structures are found between -0.1 and 
-0.2 eV along T-X and around -0.1 eV along T-M. The 



orbitally-resolved spectral functions allow to attribute 
both of these structures to the j e ff=3/2 |toj|=1/2 orbital, 
whereas the hole-like a-pocket around T is of j e ff=3/2 
\nij\=3/2 type. The two other pockets, Pm and j3x, are 
(mostly) of type j e g=l/2. 

We find that Sr2Rh04 is close to the Mott transition 
(which would occur for interaction strengths of F°=1.8- 
2.0 eV while keeping J fixed). The occupancies of the 
ieff=l/2 and j e ff=3/2 |toj|=3/2 spin-orbitals are 0.63 and 
0.89, resp. The quasi-particle weights for j s g=l/2 and 
joff=3/2 |TOj |=3/2 have been estimated as 2 , 1 / 2 =0.5 and 



J 3/2,|3/2 



=0.6. 



In conclusion, we have shown that only the coopera- 
tive effect of SOC, lattice distortions and Coulomb corre- 
lations drives the 5d compound Sr2lr04 insulating, due 
to a complete spin-orbital polarization resulting in an 



effective half-filled one-orbital (more precisely, two-spin- 
orbital) system. The isostructural and isoelectronic 4d 
compound Sr2Rh04 has smaller SOC and Coulomb in- 
teractions, leading to a less dramatic reduction of spin- 
orbital fluctuations and, as a result, a two-orbital (or four 
spin-orbital) metal. The structural and electronic simi- 
larities between Sr2lr04 and La2Cu04 - both are Mott 
insulators with one hole in an effective one-orbital sys- 
tem - together with the pronounced difference in their 
magnetic structure, suggest that the transport properties 
of doped Sr2lr04 may give valuable information about 
the role of magnetic fluctuations in high-T c superconduc- 
tivity. Oxygen-deficient Sr2lr04 is not superconducting 
[29j. Our electronic structure calculations show, how- 
ever, that doping by impurities with extremely weak 
SOC should act as giant perturbations to the spin-orbital 
structure, in analogy to introducing a magnetic impurity 
into a system with weak SOC. It would be most interest- 
ing to introduce carriers into the system while conserving 
as much as possible the very specific electronic structure, 
e.g. by doping with heavy atoms such as Os or Re. 
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